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Abstract 

We propose a definition of the running coupling constant in a SU(2) 
lattice gauge theory with twisted boundary conditions. It is based on 
the correlation of Polyakov loops extended in a twisted direction at a 
distance which is a fixed fraction of the totale lattice size. We make 
the perturbative calculation which connects this definition to standard 
regularization schemes. We find Axwisted-Polyakov/A^jg = 1.6136(2). 
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1 Introduction 



The precise determination of the strong coupHng constant ctg from expe- 
riments at intermediate energies suffers from systematic errors due to our 
ignorance of the next to leading terms of the perturbative expansion. At high 
energies, the perturbative estimates are more rehable, but even a relatively 
precise value obtained for ctg when translated into a value for A leads to a 
rather approximate result. The reason is the asymptotically free character 
of the theory which amplifies the relative error on A with respect to the 
one on ctg by a factor proportional to its inverse. Lattice calculations can 
in principle relate the determination of ctg to non perturbative quantities [1] 
with systematic errors which can be kept under control. The ideal calculation 
would be the reconstruction of the running coupling constant from a low 
energy scale normalized by a non perturbative quantity - the proton mass, 
for example - to a high energy scale where the connection to a standard 
regularization scheme like MS can be safely performed within perturbation 
theory. The highest energy scale available on the lattice is provided by the 
momentum cutoff 7r/a where a is the lattice spacing. 

However, if one wants to reconcile a sizeable total volume with a per- 
turbative scale at distances of the order of the lattice spacing he needs a 
number of lattice points exceeding the computing power available in a near 
future. A finite size renormalization group technique can solve this problem: 
this method has been widely used in the study of critical points of statistical 
systems [2] and in the context of the determination of ctg by the authors 
of references [3, 4, 5]. The coupling constant running from high to low 
energies is constructed through the following steps: 

• one defines a finite volume coupling constant renormalized at a scale 
depending upon the total volume. Its relation to the lattice bare cou- 
pling constant should be well defined and calculable in the continuum 
limit . 

• one calculates the variation of the value of the coupling between a 
volume with Ni points per side (the "small" volume) and one with 
N2 = 2Ni points (the "large" volume). This calculation is repeated for 
several lattice spacing at fixed physical volumes and extrapolated to 
the zero lattice spacing limit. 
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• in order to proceed to larger physical volumes, one needs a readjust- 
ment of the lattice spacing. This is done by matching the value of the 
renormalized coupling on a "large" volume with the one calculated on 
a "small" volume at a different value of the bare lattice coupling. 

• the last two steps are iterated until the total physical volume is large 
enough to guarantee the absence of finite size effects in the non pertur- 
bative calculation which will normalize the volume in physical units. 

In this paper we present in the case of SU(2) with twisted boundary con- 
ditions a perturbative calculation which connects a definition of ctg based on 
the correlation of Polyakov loops with the MS scheme. Section 2 is devoted 
to a short discussion of twisted boundary conditions and of the proper defi- 
nition of a twisted Polyakov loop correlation. Section 3 presents the details 
of the calculation and a comment on the absence of linear divergences in 
the observable used for defining ctg. The last section contains the final re- 
sult for the ratio Aj^jg/ ATwisted-Poiyakov where ATwisted-Poiyakov (= Atp in the 
following) is the A parameter of our scheme. 

2 Twisted boundary conditions 

The method proposed relies on a definition of ctg at a scale fixed by the total 
finite volume. The perturbative calculation which connects this definition 
to a standard one should then be performed on a finite volume. This com- 
putation cannot be done in a straightforward way because of the existence 
of gauge inequivalent field configurations of degenerate minimal energy: the 
torons [9]. They are parametrized by constant and diagonal link matrices 
which give the same contribution to the lattice action as the unity matrices 
which parametrize the standard perturbative vacuum. A solution to this 
problem is to adopt twisted boundary conditions at least in two directions 
[7, 8] which make the toron contribution to the energy higher than the vac- 
uum one. 

Normal periodic boundary conditions can be seen as the surrounding of 
the original lattice with replicas in all directions. Twisted boundary condi- 
tions define the configuration in a replica which lies in a twisted direction as 
the original one after a particular and in our case global gauge transforma- 
tion. 
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By calling Q^, the SU(7V) matrix corresponding to the twist in the u 
direction and U^{r) the link variable connecting the site r to the site r + /i 
(where r = [x,y,z,t)), a minimal twisting condition reads: 

U^{r + vL) = fi^C/^(r)fit v = x,y 

U^{r + uL) = U^{r) u = z,i (1) 

with U^{r + uL) a "replica-link", i.e. a link in a site displaced by the lattice 
size in the u direction and Qg., Qy two non commuting unitary matrices. 
Throughout the paper we will adopt lattice spacing units (a = 1). If the link 
variables U^{r) G SU(7V) are parametrized by: 

U^{r) = exp {igoA^{r)) with 

= Mr) Tr(A,(r)) = (2) 

the boundary condition (1) implies for gauge fields A^[r): 

A^{r + vL) = n^A^{r)nl i> = x,y 

A^{r + i>L) = A^{r) i> = z,i (3) 

The value of the link variable in a position which is reached after two inde- 
pendent twisted translations cannot depend upon the order in which they 
are performed. This is realized if the satisfy: 

= zn^n, (4) 

where z = exp(27rifc/7V), fc = l,...,7V — l;in our SU(2) case k = 1. 
After a combined translation the transformed link variable is: 

U^{r + xL + yL) = fi,fiyC/^(r)fiJfit (5) 

If the condition (4) is satisfied, exchanging the order of the translations 
corresponds to interchanging 0,^ and Q,y in expression (5) and leads to the 
same result. The gauge transformation for the original links 

[/» ^ K{r)U^{r)K\r + ii) (6) 

implies that the replica-links transform with a twisted gauge matrix: 

A(r + vL) = fi,A(r)fit (7) 
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In the standard Wilson action with twisted boundary conditions the pla- 
quettes at the boundary of the twisted directions contain a repHca-Hnk. For 
example, at the boundary y = L — 1, the expression of a plaquette in the x-y 
plane in terms of the original link variables is given by: 

UP' = U^{x, L-l,z, t)Uy{x + 1, L - 1, i) X 

n,Ul{x + 1, 0, z, tplUlix, L-l,z,t). (8) 

Toron configurations are no longer degenerate in energy with the usual 
vacuum because they are not transformed into themselves by the twisting 
conditions and give a contribution to the boundary plaquettes which is dif- 
ferent from the one of the standard vacuum. 

Our definition of ctg is in terms of the correlation of Polyakov loops. These 
are defined as the trace of the product of links in a fixed direction and with 
the same position in the hyperplane orthogonal to that direction. In order 
to obtain a gauge invariant operator, if the direction of the Polyakov loop is 
a twisted one, one has to trace with the corresponding twist matrix: 

L-l 

P4y,z,t) = Tr([n U,{x = j,y,z,t)]n,) (9) 

3=0 

The correlation of two twisted Polyakov loops starts at order and is a 
natural candidate for defining ctg. We take the following quantity: 

{E.,yP.{0,0,0)P,{x,y,t = L/2)) ^ ' 

where the numerator and the denominator are the ground state expectation 
values of the correlations of twisted and normal Polyakov loops respectively. 
The phase multiplying the correlation ensures the invariance under transla- 
tions in all twisted directions. 

The numerator and the denominator of equation (10) are affected by 
linear divergences which spoil their continuum limit. Various authors [10] 
have argued that these divergences do exponentiate and factorize; therefore 
the ratio of eq. (10) is a well defined quantity affected only by logarithmic 
divergences which can be reabsorbed in the renormalized coupling constant. 
In the next section we explicitely verify this statement up to order g^. 
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We define the renormalized running coupling at the scale L, g^[L), through: 

0{L) ^ T{L)g\L) (11) 
where T{^L) is the lowest tree-level contribution to O. 

3 The SU(2) perturbative calculation 

We start from the SU(2) pure gauge action: 

^ = AETr(l-C/p) (12) 
yo p 

where we remind that when the links are in the replica they are replica-links 
satisfying the twisted boundary conditions (1). 

The connection of the definition in eq. (11) to other regularization schemes 
requires the calculation of the one-loop correction to the observable O, i.e. 
of terms up to order g^. These come from the expansion of the action up 
to order g^ and of the operator up to order g^. We have performed our 
calculation in configuration space to which refer our Feynman rules. The 
explicit form of the various terms defined by: 



S = S,+ 5o'^meas + 50^^^ + qIS'^^ + ^O^hLt + ^O^host + ^{gl) (13) 

is given in Appendix A. The Feynman diagrams are represented in figu- 
re 1 and 2. 

The basic building block is the gluon propagator which, due to the twisted 
boundary conditions, has peculiar properties which depend upon the explicit 
form chosen for the Q, matrices. For our SU(2) case, the non trivial fi's can 
be identified with the Pauli matrices: 

= io-i = io-2 (14) 

If the gauge fields of the SU(2) algebra are expressed in terms of the Pauli 
matrices cr^ (a=l,2,3): 

Mr) = E Alir) (15) 

a=l 
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equation (3) reads: 



Alir + vL) = Alir) 



u = x,y 



u = z,t 



(16) 



which, together with eq. (14), implies that the gauge fields of colour 1 are 
periodic in the direction y over an anomalously enlarged space of size 2 X L. 
The same happens for those of colour 2 in the direction x and those of 
colour 3 in both x and y directions. The propagator is diagonal in colour 
and shares the periodicity properties discussed above as a function of the 
distance between two points. 

Equation (14) associates to each twist matrix a corresponding colour: in 
our basis the one corresponding to translations in direction x (y) is associated 
to colour 1(2). One can define new vector potentials which have a standard 
periodicity in all directions. They are related to the old ones by a local 
colour-depending phase, according to: 



ni = (0,1,0,0) 712 = (1,0,0,0) 713 = (1,1, 0,0) (17) 



Operators built from links expressed in terms of the new vector potentials 
have normal translation invariance properties. The phase appearing in the 
definition of the correlation between twisted Polyakov loops is just the net 
result of the transformation from a loop expressed in terms of the new links to 
one in terms of the old. It can also be verified explicitely that tracing a single 
twisted loop with the appropriate twist matrix times a phase corresponding 
to the associated colour produces an ordinary translation invariant operator. 

The calculation in configuration space was done using FORM [12], an 
algebraic manipulator, in order to deal with the plethora of Wick contrac- 
tions. The sums over the positions were done with a numerical FORTRAN 
program. Their number in the most time consuming diagram was reduced 
from twelve to five. The gluon and ghost propagators in configuration space 
were obtained by doing a fast Fourier transform of their analytic expression 
in momentum space (see Appendix B for details on the momentum space 
formulation). 



Alir) = Alir) exp(-i-(7i„ • r)) 



with 
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4 The results 

Our calculation passed successfully some crucial tests. 

(i) The diagrams contributing to the gluon self energy exhibit ultraviolet 
quadratic divergences which appear as contributions growing quadra- 
tically with the number of lattice points per side. These divergences 
cancel. 

(ii) The sum of diagrams of the twisted Polyakov loops correlation has 
genuine linear divergences which cancel against the product of the de- 
nominator at order times the tree level numerator. 

(iii) After the above cancellations one is left with the logarithmic diver- 
gences which make the renormalized ctg running with the lattice size. 
Their coefficient Sq is known from continuum calculations and is given 
by: 



127r2 

to be compared with our result 



2/3o = s'^ ~ 0.092878 (18) 



50 = 0.092880 ± 0.000005 (19) 

In the following we shall adopt this notation: 

0{L) = T{L)[gl + gML)] (20) 

where D[L) is the (tree-level normalized) one-loop contribution. 

The L dependence of the one-loop contribution can be parametrized as^ : 

D{L) =ro + so log(L) + ^ + ^ + ^ log(L) + . . . (21) 

term of the form ^log(i) is ruled out: the lattice pure gauge action differs from 
the continuum one only by terms of order , i.e. in our units of order which implies 

that terms of the above form cannot be generated at one loop level. The ri/L term could 
be produced by combining linear divergences with l/L^ corrections. However, it turns out 
to be absent because the linear divergences appearing in the numerator and denominator 
of eq. (10), which defines O, fully factorize and exactly cancel in the ratio. 
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ij-range 




So 


ri 


4— oz 


u.zoy 1 


U.iUiz 


U.io4y 


8-32 


0.3160 


9.430 • 10-2 


3.061 • 10-2 


12-32 


0.3194 


9.345 • 10-2 


1.388 • 10-2 


16-32 


0.3206 


9.316 • 10-2 


7.325 • 10-2 


18-32 


0.3211 


9.304 • 10-2 


4.381 • 10-2 


20-32 


0.3213 


9.300 • 10-2 


3.476 • 10-2 


22-32 


0.3219 


9.286 • 10-2 


-5.167-10-4 



Table 1: Results of the first fit showing numerical evidence for ri = 

Various fits of D[L) were made, for a series of lattices ranging from L = 4 
to L = 32. 

The first was addressed to obtain a numerical evidence for the absence of 
the 1/ L terms; in table 1 we show the results of the fit: 

D{L)=ro + solog{L) + j- (22) 

for different ranges of lattice sizes. While the fitted values of the coefficients 
of the logarithm Sq and of the constant ro show no appreciable dependence 
upon the smallest lattice size L used in the fit, the value of the linear term 
ri decreases systematically, providing evidence for the absence of such a 
contribution. 

Assuming that ri = 0, the second set of fits (see table 2) was obtained by 
fitting D[L) with the form ro + So\og[L) and led to a precise determination 
of the coefficient of the logarithm with the result in eq. (19). 

The precision (~ 5 X 10-^) of this determination of ro is better than what 
we expected: neglecting in the expression for D[L) terms of the form and 
\og[L) induces an error, which, in principle, should be of order 10-2, for 
the range of lattice sizes of interest. In order to clarify this point, we made (for 
various ranges of lattice sizes) various fits of the quantity D(L) — s^^log(L) 
with the expression: 

-o + 9 + ^log(i). (23) 

In the range of lattice sizes 22 < L < 32, which leads to a stable determina- 
tion of ro and Sq in the previous set of fits (see table 2), we found for r2 and 
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L— range 




•So 


C 


16-32 


0.32199 


9.2833 • 10-2 


2.8 • 10-^ 


18-32 


0.32191 


9.2857 • 10-2 


1.2 • 10-^ 


20-32 


0.32188 


9.2868 • 10-2 


1.0 • 10-^ 


22-32 


0.32184 


9.2878 • 10-2 


9.2 • 10-^ 


24-32 


0.32183 


9.2883 • 10-2 


9.0 • 10-^ 


26-32 


0.32184 


9.2878 • 10-2 


9.0 • 10-^ 



Table 2: Results of the second fit determining the value of Sq, where ^max = 
Max{i} {\D{L) - so\og{L) - ro\/\D{L)\) 

S2 values that make the combination [r2 + S2log[L)) / always smaller than 
10-^, with the value for ro stable within a precision of 10-^. In this range of 
lattice sizes we are allowed to neglect in our fits the terms of order and 
L~^\og[L) and we obtain results with a precision of some units in 10-^. 

The last set of fits was made with the coefficient of the purely logarith- 
mic term fixed to its theoretical value, s^^, to get the best determination of 
the constant term which fixes the renormalization scheme (see table 3 and 
figure 3). For various ranges of lattice sizes, we made a fit of the quantity 

D{L) - s'^\og{L) (24) 

with a single constant parameter {vq). 
Our final result is: 



ro = 0.32185(1) (25) 

where the error was estimated by emaxfo — 10-^ (cmax is defined in table 3). 

This constant settles the relation to any other regularization scheme and 
allows a comparison of the numerical result for the running coupling constant 
with other definitions. The ratio Axp/A^jg gives the relationship between our 
definition of ctg and the one of the MS scheme. 

We make use of the relation [11] between the coupling constant in the 
MS regularization scheme, and the lattice renormalized coupling, ^f^tt, 
defined by 

2/5oiog(/x/Aiattj 
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ij-range 


ro 


^max 


i.Z-6Z 


0.32188 


AO 1 r\— 4 
4.0 • 10 


io-oz 


U.ozioo 


i.U • iU 


18-32 


0.32185 


3.6 • 10"^ 


20-32 


0.32185 


1.9 • 10"^ 


22-32 


0.32185 


1.9 • 10"^ 


24-32 


0.32185 


1.9 • 10"^ 


26-32 


0.32185 


1.8 • 10"^ 


28-32 


0.32185 


1.6 • 10"^ 



Table 3: Final fit for ro, where Cmax = Max{i} [\D[L) — s^log^L) — ro|/|ro|) 

where fl^fattC-'^) = 9o- knowledge of ro fixes the relation between ^^^p and 
fl^fatt) which reads, to one-loop order: 

9tA(^-') ^ w\a ^ = + ^o^iatt(/^"') (27) 

2/3olog(/x/ATp) 

where is a physical length. The desired relationship is given by: 

^ = TZ:-P(5-) = 1-6136(2) ■ (28) 

Appendix A 

We have made the perturbative calculation of the observable O in eq.(lO) up 
to the order for the expectation value in the numerator and g^ for the one 
in the denominator. The expectation values have the standard definition in 
terms of ratios of functional integrals: 

After the Feynman gauge fixing, the functional integral becomes: 

J DcDcDAe-^^'^'''''^C[A] 
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where c, c are the ghost fields and S is given by 

S = S, + glSmeas + 5o5&^ + gls'^^ + +g.s'^Zst + alsfL, + ^(^o') (31) 

The Feynman diagrams needed to perform the calculation are given for the 
numerator of eq.(lO) in figure 1 and for denominator in figure 2. 
Diagrams A, B and C need the expansion of the twisted Polyakov loop up to 
order g^. With the exception of diagram I which is zero, the remaining ones 
all need the twisted Polyakov loop only to order g^. Moreover, diagrams D 
and E need the gauge field action up to order g^, diagrams F contains the 
measure at order g^, and diagrams G and H involve the ghost action up to 
order g^. Diagrams in figure 2 involve only the Polyakov loop at order g^. 

The expansion of the operator is straightforward. The one of the gauge 
action is a little bit more involved and we report in the following the explicit 
expressions for the quantities of eq.(31). 



S'£' =-4E.E^>{ i»xi> + i>).[X(r + /i) + X(r)]} (32) 

S^w^ =-|E.E^>{ +4 [X(r) • i> + !>)] [X(r).X(r)] + 

+4 [X(r-) • A^{r + v)] [A^{r + v) ■ A^{r + v)] + 

-4 [X(r) • X(r + fi)] [X(r) • X(r)] + 

-4 [X(r) • X(r + /i)] [X(r + /i) • X(r + /i)] + 

-2 [X(r) • A^{r)] [X(r) • X(r)] + 

-2 [X(r) . A,{r)] [A,{r) ■ A,{r)] + 

-2 [A^{r + (i) ■ A^{r + v)] [A^{r + v) ■ A^{r + v)] + 

-2 [A^{r + (i) ■ A^{r + v)] [A^{r + (i) ■ A^{r + (i)] + 

+6 [X(r) • X(r)] [A^{r + i>) • i^(r + v)] + 

+6 [X(r) • X(r)] [X(r + /i) • X(r + /i)] + 

+3 [X(r) • X(r)] [i^(r) • A^{r)] + 

+3 [l^(r + v) ■ A^{r + v)] [A^{r + /i) • A^{r + /i)] + 

+ [X(r).X(r)] [X(r).X(r)] + 

+ [X(r + /i) • X(r + /i)] [X(r + /i) • X(r + /i)] + 
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+ 12[A^{r) ■ A^{r + i>)] [A^{r + fi) ■ A^{r)] + 
+ 12[X(r) • A^{r)] [A^{r + u) ■ X(r + /i)] + 
-12[X(r) • X(r + (i)] [A^{r + i>) • A^{r)] + 
+ 12[X(r + (i) ■ A^{r)] [X(r) • Mr)] + 
+ 12[X(r + /i) • A^{r)] [A^{r + i>) • A^{r + i>)] + 
-12[i^(r + i>) • X(r- + A)] • X(r)] + 

-12[A^{r + u) ■ A^{r)] [X(r) • X(r)] + 
-12[X(r) • X(r + /i)] [i^(r + i>) • i^(r + i>)] + 
-12[X(r).i»] [i^(r + i>).i> + !>)]} (33) 

5™ea. =|E.E^ (34) 

^ilt = 4 E. E,{ +?(r-) • Air) X (c(r) + c(r + fi))] 

-g(r) • [i^(r - A) X (c(r - /i) + c(r))]} (35) 

SSost =-|E.E,{ +[?(r).(c(r)-c(r + A))][i.(r-)-i»] 

-[g(r) . A,{r)][A,{r) ■ (c(r) - c(r + /i))] 
-[g(r) • (c> -fc)- ^{r))]A{r - fc) ■ A^{r - fi)] 
+ [g(r) . A,{r - f^)][A,{r - j,) ■ (c> - /i) - c»)]X36) 

Vectors are defined in the colour space and scalar and vector products 
accordingly: 

A = (AW,A(^),A(^)) 
c = {c'<^\c^^\c^^^) 

The sums Er range over all the points of the original lattice. This implies 
that some terms contain replica vector fields which are related to the original 
ones by eq.(16). 
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Appendix B 



In this appendix we want to discuss the relation between the colour degrees 
of freedom and the extra momentum degrees of freedom introduced on the 
lattice by the twisted boundary conditions (see equations (1) and (3)). 
Expanding the SU(7V) gauge field A^[r) into plane waves 

^» = ^Er.i.We^'='-+^^ (37) 
^ k 

the twisting condition (3) implies: 

n^T^ni = e^'^-'^Tk u = x,y (38) 

where Fj. is a complex N X N matrix. It can be shown (see ref. [6]) that 
a non-zero solution of equation (38) exists if and only if the momentum 
components in the twisted directions [u = x,y) satisfy 

= + k^ with 

^^ = ^< < = 0,l,...,Ar-l (39) 

These momentum components, as expected, are quantized, but the k±^ term 
in the momentum adds to the ordinary momentum degrees of freedom those 
of a lattice N times larger in each twisted direction. In non-twisted directions 
(i) = i,i) one can think of k^, as a sum k^ + fcj", with k^ = 0, because in 
these directions the lattice of size L exhibits the ordinary periodicity. 

Up to an overall arbitrary phase, the solution of equation (38) is unique 
and reads 

Tk = fi;<fi^^^z^("^^+<+^)("^^+<) (40) 
where we used the Q property 

n^ = {-lf-'l = x,y. (41) 

From eq.(40) follows that the Fj-'s are SU(7V) matrices, that they depend 
only upon k± and that [5] their normalization is: 
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^Tr(rt,r,.) = 4V (42) 

where is 1 if, for i/ = x,y, = n^(mod N) and zero otherwise. 

Equation (37) can be rewritten as 

Mr) = j-, E ^,.i,(F^fc^)e«^ (43) 
where the condition Tr(yl^(r)) = is implemented by requiring 

i^(F^,fc^) = if n^ = 0(mod7V) . (44) 



The gluon propagator, in momentum space and in the Feynman gauge 

is: 
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where k = k^^^ + k-^, 

Xk- = l-6,.,o , k' = Aj2sin'{^) (46) 
and the bilinear [k-^,k-^) is defined as 

{k^ , k^) = na^Ua: + fiyUy + {fia: + ny){na: + Uy) . (47) 

The comparison of equation (43) with the ordinary Fourier and colour 
decomposition of the gauge field A^[r), shows that the colour degrees of 
freedom of Af^[r) are transferred into the k-^ (unphysical momentum) degrees 
of freedom: their number is just — 1, due to the trace condition (44). 

We now specialize to our particular case, SU(2), with the twist matrices 
Qg., Qy given by equation (14). 

^In this gauge, the ghost propagator in momentum space has the same form as the 
gluon one, with replaced by 1. 
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It follows that: 

rfc|(nJ-=0,nJ-=l) = -^O-l 
rfc|(nJ- = l,nJ-=0) = +^0-2 
rfc|(n^ = l,n^ = l) = -io-3. (48) 

By rewriting equation (43) in the form 

Mr) = j-, E ^,.i,(F^fc^)e«^ = ^E^ci^(^"V"■+'^, (49) 
we find: 

i;,(0 = -ii,(F^fc^)l(„.=o,n,-=l) 

i^(F^) = +ii^(F^fc^)|(„.=l,„.=o) 

i^(F^) = -ii^(F^fc^)|(„.=l,„.=l), (50) 

where a correspondence between the k-^ and the colour c indices is explicitly 
established. The details of this correspondence depend on the choice for the 
Q matrices. 

The gluon propagator (45) becomes in this colour basis and for the SU(2) 
case : 

An ordinary Fourier transform gives the corresponding propagator in con- 
figuration space (which is diagonal in colour indices): 

G%{r) = ^ E G'^,(F^)e^'=-^'-+^'="(^)'- (52) 

where k-^{c) is given by the correspondence established in equation (50). The 
phase factor containing k-^[c) is responsible for the peculiar, colour dependent 
properties of G^^(t") under translations by L in the twisted directions. 
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